BI-TP 2012/06, ULB-TH/12-06 



Hairy Black Holes in Massive Gravity: Thermodynamics and Phase Structure 

Fabio Capela 

Service de Physique Theorique, Umversite Ltbre de Bruxelles (ULB), 

CP225 Boulevard du Triomphe, B-1050 Bruxelles, Belgium 

and Germano Nardini 

Fakultdt fiir Physik, Universitdt Bielefeld, 

Universitatsstrafle 25, D-33615 Bielefeld, Germany 



o 

< 



The thermodynamic properties of a static and sphericaUy symmetric hairy black hole solution 
arising in massive gravity with spontaneous Lorentz breaking are investigated. The analysis is car- 
ried out by enclosing the black hole in a spherical cavity whose surface is maintained at a fixed 
temperature T. It turns out that the ensemble is well-defined only if the "hair" parameter Q char- 
acterizing the solution is conserved. Under this condition we compute some relevant thermodynamic 
quantities, such as the thermal energy and entropy, and we study the stability and phase structure 
of the ensemble. In particular, for negative values of the hair parameter, the phase structure is 
isomorphic to the one of Reissner-Nordstrom black holes in the canonical ensemble. Moreover, the 
phase-diagram in the plan (Q, T) has a line of first-order phase transition that at a critical value 
of Q terminates in a second-order phase transition. Below this line the dominant phase consists of 
small, cold black holes that are long-lived and may thus contribute much more to the energy density 
of the Universe than what is observationally allowed for radiating black holes. 
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I. INTRODUCTION 

The discovery of the accelerated expansion of the Uni- 
verse [T] has revived interest in the theories of gravity 
that can explain such observation without invoking dark 
energy. One class of these models is called "massive grav- 
ity" . Models of massive gravity fulfilling Lorentz invari- 
ance have been recently constructed [2] and proved to be 
free from ghosts and instabilities at full non-perturbative 
level Oil]. However, there are very few known solutions 
for such models. Instead, models of massive gravity with 
Lorentz violation are better understood at present. In 
these scenarios the spontaneous breaking of Lorentz sym- 
metry is achieved by scalar fields that are coupled to grav- 
ity in a covariant way through derivative couplings [5HZ] ■ 
As a consequence of this interaction, Lorentz violation is 
transmitted to the gravitational sector and the graviton 
becomes massive j6- 8' . In a wide region of the parameter 
space these Lorentz breaking models are, at perturbative 
level, free from ghosts and tachyonic instabilities around 
Minkowski |H] and curved backgrounds [5] . They can also 
exhibit infrared modifications of the gravitational behav- 
ior [£,. Moreover, they are expected to reach the strong 
coupling regime at energies that are much higher than in 
Lorentz-invariant models of massive gravity |6H9], even 
though it has not yet been rigorously proved. Interest- 
ingly, these Lorentz breaking theories are formulated in 
a non-perturbative way, making the study of black holes 
possible [TOHTSj. On top of that, black hole solutions are 
far richer than in General Relativity (GR) due to the 
presence of "hair" parameters [TTl [T^ . 

In GR the existence of black hole solutions raises ques- 
tions of consistency with the general laws of thermody- 
namics. For instance, it has been argued [I^ that black 
holes have temperature and entropy. The same conclu- 



sion was apparently reached by means of path integral 
methods P^. However, the proof turned out not to 
be fully consistent because of a technical subtlety: the 
path integral approach requires black holes in thermal 
equilibrium with its surroundings and this situation is 
never fulfilled by Schwarzschild solutions. Instead, anti- 
de-Sitter (AdS) black holes achieve thermal equilibrium 
and the path integral procedure proposed in Ref. [15] 
can be consistently applied HHj. Moreover, it turns out 
that black holes are just one of the phases permitted in 
AdS space and first order transitions between a black 
hole and globally-AdS spacetime may occurjj On the 
other hand, it might be expected that a similar result 
also arises for Schwarzschild black holes if one mimics 
the AdS cosmological constant by an infrared cutoff. In 
fact, as it was proven by York 19], after having enclosed 
the Schwarzschild black hole inside a spatial spherical 
surface at a fixed temperature, the configuration reaches 
thermal equilibrium and its phase structure turns similar 
to the AdS one. 

It is interesting to understand whether one can apply 
the York's approach to investigate thermodynamics of 
hairy black hole solutions in Lorentz breaking massive 
gravity. This issue is the subject of the present paper and, 
as we will see, the conclusion depends on whether the 
black hole "hair" parameter Q is conserved in the cavity. 
In particular, when Q is assumed constant rj we obtain 
a rich thermodynamical phase structure that includes. 



^ In the context of the AdS/CFT correspondence 1171 . this transi- 
tion turns out to be dual to the confinement-deconfinement phase 
transition in large- A^ gauge theories [18j. 

^ In practice, we expect that the results obtained under this hy- 
pothesis can be extended to situations where Q varies very slowly. 



for instance, first and second order phase transitions and 
phases consisting of stable and small, cold black holes 
(plausible dark matter candidates). 

The paper is organized as follows. In Sect. Ill] we briefly 
review the model of massive gravity and we sketch the 
eventual difficulties that may arise when applying the 
York's procedure. In Sect. |III| we review the static, spher- 
ically symmetric and asymptotically-flat black hole solu- 
tion of massive gravity and we check whether it can be 
embedded in a spacetime that has a periodic Euclidean 
time coordinate and is bounded by a spatial sphere. This 
compatibility requirement is needed in order to apply the 
York's procedure, which we carry out in Sect. IV assum- 
ing Q to be conserved in the ensemble. In particular, we 
first regularize the black hole action by the "subtraction 
background" method [TS] , then we enclose the black hole 
solution in a sphere at a fixed temperature and, finally, 
we obtain the on-shell action of the ensemble. As we 
prove in Sect. W\ this procedure provides a well-defined 
ensemble that we use to investigate black hole thermody- 
namics. It turns out that, depending on the value of Q 
inside the cavity, black holes can evolve in qualitatively 
different ways. In any case, their evolution always re- 
spects the Bekenstein-Hawking formula. In Sect. |Vl| we 
relax the assumption concerning the conservation of Q 
and we see that in such a case the York's procedure is 
inconsistent. Finally, Sect. |VII| is dedicated to summa- 
rize the main results of the paper and Appendix contains 
some analytic expressions valid for a specific choice of 
massive gravity parameters. 



II. MASSIVE GRAVITY AND THE 
SADDLE-POINT APPROXIMATION 

Massive gravity is conventionally described by the ac- 
tion IMHl 
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dl^Xy/g 
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where Latin (Greek) indices run on space (spacetime) 
components. The first integral in Eq. (IT]) is evaluated 
on the manifold M. with metric g and contains two con- 
tributions: the usual Einstein-Hilbert term and a func- 
tion T of four scalar fields (/)'' that are minimally cou- 
pled to gravity by covariant derivatives. The second in- 
tegral is instead the Gibbons-Hawking- York boundary 
term [15] f^Oj, where 7ah is the metric induced on the 
boundary 9A^ and K is the trace of the extrinsic curva- 
ture Kij — ^^^VkTij of dM with unit normal n*. Such 



a boundary term is required to have a well-defined vari- 
ational principle in the presence of the border (9A^ . 

The action (II]) describes a low-energy effective theory 
valid below the ultraviolet cutoff A, which perturbative 
analyses estimate to be 0{\/rnMp\) [SHS], where Afpi = 
•\/l/87r and m are the Planck and the graviton masses, 
respectively. Its vacuum flat-spacetime solution has the 
form 



Qiif — VtJ.v 



C = A^i 



"flat 



= A^ 



(2) 



which induces a spontaneous breaking of Lorentz symme- 
try. The background (J2]) preserves rotational symmetry 
when the function T is invariant under rotations in the 
internal space of the fields cj)^. Moreover, the action (IT]) 
is invariant under the symmetry 



^• 



e' 



(3) 



where 0' are arbitrary functions of (jP . This symmetry 
ensures that perturbations around the vacuum contain 
only two propagating degrees of freedom ^ , correspond- 
ing to the two polarizations of a massive graviton. 

To study the thermodynamics of this model, we use 
the Euclidean path integral 



Z = JvgVcj) exp(-/[5,,^]) , 



(4) 



which is evaluated by integrating over all metrics and 
scalar fields satisfying particular boundary conditions. In 
the semi-classical limit Z is dominated by the stationary 
points of the action. This can be checked by expanding 
the path integral around a classical solution. Indeed, if 
the expansion provides a leading term that is finite, a 
linear term that vanishes on-shell, and a quadratic term 
that is positive definite, then the function Z can be ex- 
pressed as 



'-f[ffcl,0cl] 
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and can be interpreted as the partition function of the 
model. However, such a derivation is not straightforward 
for the action (IT]) since the three properties listed above 
might not be fulfilled. In fact: 

1. The on-shell leading term of the action diverges. 
This is a familiar problem in general relativity that 
is addressed using the "background subtraction" 
regularization technique [15j . Following this sub- 
traction scheme, we take the vacuum solution ^ 
as background {go, (/iq) and we consider Ie as reg- 
ularized action, defined as 



Isig, (f>) = l{g, 0) - /(go, ^o) 



(6) 



In this way Ie is finite for the class of fields 
{g, (p) that asymptotically approach the background 

(50 >o); 



2. The linear term may eventually not vanish for all 
perturbations around the classical solution. The 
non-vanishing behavior of such a term comes from 
the boundary contributions: 
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For the action (fTl) tt"'' and ttJJ are given by 
^ -[K 

dx dW'i x^ f" ^ 
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where ria is the outward pointing unit normal to 
dM. Therefore, some boundary conditions on the 
fields g^i/ and (/)'^ should be imposed in order to have 
a well-posed variational problem, i.e. 5I\^^ = 0; 

3. The Gaussian integral in Eq. (Isl) corresponds to the 
one-loop contribution. Such contribution may di- 
verge since the integration involves the determinant 
of second-order elliptic operators that cannot be 
regularized when negative eigenvalues are present. 
Indeed, if there exist negative eigenvalues, the den- 
sity of states grows so rapidly that the ensemble 
turns out to be ill-defined and Z cannot be used 
to determine thermodynamical quantities. For this 
reason, in order to be able to study the thermody- 
namical properties of black holes, we need first to 
stabilize the ensemble. This can be performed as in 
Ref. [19]: we place the black hole inside a surface 
maintained at a fixed temperature. 

In the next sections we explicitly show how to imple- 
ment these three procedures for black hole solutions in 
massive gravity. 



III. THE BLACK HOLE SOLUTIONS 

By the coordinate transformations r' -^ r ^ r{r') and 
t' ^f t ^ t-l-r(r'), the generic ansatz for the static spher- 
ically symmetric solution in Euclidean spacetime can be 
written as [TH HI] 

ds^ = a{r)dt^ + p{r)dr'^ + r^ [dO^ -|- sin^ dd^p^) , 
0° = A2H^ + %)] , (10) 

0^ = 0(r)ll^. 
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For asymptotically-fiat black holes the solution of Ein- 



stein and Goldstone's equations, respectively E^ 



(1) 







and yP' = 0, has a known analytical expression if the 
function F takes the form [T^ 

12 /2 \ 

F= T^ + 6[ T + l) wi-iA'i+SwiWa- 2^3 + 12, (13) 

where A is a positive constant and Wn = Tr(Ty"). In such 
a case, for A ^ 1 , the ansatz ( 10 1 provides the black hole 
solution 

a[r) = 1 ^ , 

p{r) = 4t : (14) 



a{r) 



h{r) 



± 



dr 
a 



Q A(A-l) , ^ 

1 - a — — ^ rr^ h 1 



12m2 



rA + 2 



1/2 



(j){r) = r , 

which depends on the two arbitrary integration constants 
M and Q. In the following we restrict our analysis to the 
class of solutions (14 1 with A > 1, so that the gravi- 



tational potential is asymptotically Newtonian and the 
parameter M coincides with the Arnowitt-Deser-Misner 
(ADM) massR Moreover, we forbid naked singularities, 
i.e. a{r) must have real roots and the largest of them 
determines the radius of the event horizon. Depending 
on the signs and relative values of the parameter M and 
Q, the following cases arise: 

9 Q > and M > : The existence of the hori- 
zon is guaranteed: at all distances the gravitational 
potential is attractive and stronger than (or, for 
Q = 0, equal to) the usual Schwarzschild black hole 
potential. Therefore r^ is never smaller than in the 
standard case. 



For black hole solutions in massive gravity the explicit 
expression of this ansatz can be obtained by imposing 
Eqs. ( 10 1 to fulfill the equations of motions of the ac- 
tion (llj). The black hole solution has then to be an ex- 
tremum of the variation of the action with respect to the 



^ These features have been discussed in Ref. 12 for A < 1. For the 
particular case A = 1, yielding a different solution from | |14[ |, see 
Appendix F of Ref. Ill| . Such a solution produces a divergent 
ADM mass . 



Q > and A/ < : The Newton's potential 
is repulsive at large distances and attractive near 
the horizon. This possibility is interesting because, 
even for Q = 0, it does not have a corresponding 
case in GRq Nonetheless, we do not analyze such 
negative mass configurations since they are incom- 
patible with black hole thermodynamical laws [13] . 
As stated in Ref . [23 , it seems likely that the cause 
of these incompatibilities is not Lorentz violation 
but the existence of negative energy states. 



• Q < and A/ > 
the condition 



The horizon only exists when 



2Af > A|Q| 



l/A 



1 

A-1 



(15) 



is fulfilled. In this case the Newton's potential is 
always attractive until reaching the horizon but the 
attraction is weaker (which makes the event horizon 
radius r+ smaller) than in the Schwarzschild case. 



A. The boundary conditions 

As it has already been mentioned in Sect. |llj the pro- 
cedure to study equilibrium thermodynamics requires to 
enclose the asymptotically-flat black hole within a finite 
volume surface and then send the surface to infinity |19| 
(see Sect. |lV]for details). In the following, we consider 
a spherical cavity of radius r^ as the boundary and, to 
analyze the system at finite temperature, we impose pe- 
riodicity on the Euclidean time. Due to the presence 



of this border, the black hole solution (14) must fulfill 
some boundary conditions in order to satisfy (5/|ci = 0. 
The momentum ttq conjugate to the scalar field 0" on 
dAi vanishes when evaluated at the solution (14). Then, 
there is no need to fix the scalar field 0" at the bound- 
ary. Analogously, as a remnant of the internal spherical 
symmetry of the scalar fields </)', i.e. {<j)'^)^\dM = ''"B, 
no boundary conditions emerge on the scalar sector. In- 
stead, the momenta tt"'' conjugated to the induced metric 
Jab are not null when evaluated at the background solu- 
tion ( 14 ) . For this reason we conclude that we only need 
to impose 6jab = to have vanishing boundary terms. 
Further conditions arise due to the compactification of 



* The choice M < is problematic in GR where: i) it leads 
to naked singularities; ii) it violates the null energy condition 
Tfiiyk'^k'^ > (being fc^* a future-pointing null vector field), which 
holds for the matter stress tensor |21| and implies the positivity 
of the ADM mass. Neither of the arguments exist in massive 
gravity: i) at short distance the repulsion changes to attraction, 
which creates the event horizon; ii) the stress-energy tensor of 
the scalar fields T^j^ does not satisfy the null energy condition, 
allowing for negative mass states to be constructed, e.g., as in 
the ghost condensate model 1221 . 




FIG. 1. The Hawking temperature Th as a function of the 
event horizon radius r+ for positive and negative scalar charge 
Q and for the conventional Schwarzschild case Q = 0. For 
Q < the Hawking temperature reaches a maximum at r+ = 

\2 iM^ni/A 



r-„,ax= [i\''~l)\Q\] 



the Euclidean time. Regularity of the metric at the hori- 
zon fixes the periodicity a,s t ^ t + Ph , where /3h is 
related to the Hawking temperature Th by the equality 
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Ph 



dra 



inrA 



r+ 



1 + (A-1)- 



(16) 



Interestingly, the Hawking temperature behaves differ- 
ently for positive and negative scalar charges (see Fig. [I]). 
For Q > the temperature is larger than for Q = and 
decreases with the increasing of the event horizon radius 
r+ , as in the case of the conventional Schwarzschild black 
hole. For Q < the temperature is smaller than in the 
Schwarzschild case and its non-monotonic dependence on 

r_|_ reaches a maximum at r+ — rmax = [(A^ — 1)|Q|] 
The periodicity of t also imposes a periodicity on the 
scalar field ((^'{t^r) ~ (jP{t + /3_f/,r) which is linearly de- 
pendent on time. This does not provide any further con- 
straint on the solution ( |l4| ) since the periodicity of (jP is 
allowed due to the presence of the global shift symmetry 
(jp ^ (jp + C in the action ^. 

Notice that the Hawking temperature measured by an 
observer at the position r is given by 



T{r) 
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1 + (A 



1)^ 
' + 



Attta 



(17) 
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In case of thermal equilibrium, setting the temperature 
measured at the surface r — r^ univocally determines the 
temperature of the configuration contained in the cavity. 
Once one has set rg and T{rig) the boundary dM is 
unambiguously fixed. 



IV. THE ON-SHELL ACTION 

In this section we briefly review the background sub- 
traction procedure needed to obtain a finite on-shell ac- 



tion. Subsequently, we compute the Euclidean action of 
the black hole solutions ( 14 ) in terms of the boundary 
conditions (T(?'b),?'b) under the assumption that Q is a 
conserved quantity in the cavity. 



A. The regularization procedure 

A priori, the Euclidean action (II]) may present an in- 
tegration problem at r = r_|_ and another at r — > oo. The 
former is solved as usual by assigning the period j3fj to i, 
as in such a case the metric ( 14 1 extends smoothly onto 
the event horizon. The latter is instead more cumber- 
some. To regularize it 19 , we perform the integration 
from r = r_|_ up to the infrared cutoff r = rg, we subtract 
off the action of the vacuum flat space (plF] and finally 
we send rg to infinity. In this way the Euclidean action 
is regularized and defined as in Eq. ([6|. 



with 



Notice that in a spherical cavity of finite volume the 
on-shell action (with the metric regularized at r = r_|_) 
is always finite and, in addition, the thermodynamic sta- 
bility is guaranteed [19]. Hence, as long rg is finite, in 
principle there is no need to subtract any background to 
make the action finite. Instead, this procedure is nec- 
essary to have no divergent thermodynamical properties 
when the limit rg — >■ oo is taken. Then, in view of this 
limit, we apply the regularization procedure as a first 
step to analyze black hole thermodynamics. 



B. The regularized action for hairy black holes 



A = Anrl , 



Sby = rs 



1- v^a(re) 
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8^ 



(20) 



d^x^/a [k — ko) 



Bi 



where k is the trace of the extrinsic curvature of the two- 
boundary Bi = dA4\t=i , (7 is the induced metric on Bi 
and ko refers to the extrinsic curvature of Bi embedded 
in the vacuum space Q. 



The last equality in Eqs. ( 20 ) shows explicitly that £by 
is the Brown- York quasi-local energy [21 



corresponding 

to the Hamiltonian that generates the time translation 
at the two-boundary Bi. For this reason we can interpret 
£by in Eq- ( 19 ) as the energy of the black holer] On the 
other hand, plugging Eq. (14) in Eq. (20) reads 



M = £ 



BY 



'--BY 



2rB 2r 



x-i 



(21) 



showing that the ADM mass M is the total energy of the 
black hole in the limit rg — > cx). For this reason, when we 
send the cavity surface to infinity, the regularized on-shell 
action is given by 



Ie = I3hM - -Abh 



(22) 



which still vanishes when the event horizon goes to zero as 
it reproduces the action of the subtracted vacuum space- 
timeQ A detailed discussion about the finite energy of 
the black hole solutions (14 1 can be found in Ref. [25] . 



An apparent ambiguity may arise in the procedure 
to regularize the action (IT]): since the subtracted back- 
ground (J2]) is regular everywhere, it does not require 
any specific periodicity of the time coordinate. However, 
black hole and background metrics have to match at the 
boundary surface r — r^. Thus, the time periodicity /? 
of the background has to be 



/? = T{rB)- 



(18) 



where T{rBy^ ^ l^nVairB) as in Eq. pj] ). 

Once the action (JTI has been regularized, the evalu- 
ation of Ie for the black hole solution (JmI with radial 
integration r+ < r < rg is straightforward. Indeed, as 



the metric ( 14 ) and J- of Eq. ( 13 ) are stationary, the time 



integration in Ie just gives rise to a multiplicative fac- 
tor /3. Moreover, the integration on the other variables 
can be performed in a closed manner. The regularized 
(quasi-local) action of the black hole then results to be 
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(19) 



V. THERMODYNAMICS AND PHASE 
TRANSITIONS 



Eq. (17) evaluated at r = rg produces a function of 



the ADM mass M in terms of the parameters of the en- 
semble /3, rg and Q. Depending on the particular values 
of these parameters, there can exist zero, one or multiple 
black hole solutions that are allowed inside the cavity. 
When several configurations are permitted, phase tran- 
sitions may occur. In the present section we analyze this 
issue under the hypothesis that the scalar charge is con- 
served inside the cavity (for considerations without this 
assumption see Sect. VI). 



A. The phases 

In order to determine the number of black hole solu- 
tions allowed inside a cavity containing a given scalar 



^ Of course, this subtracted background is not I 
tlie periodicity of time is maintained. 



Remind that by construction Eby = for the vacuum spacetime 
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^ Of course, this conclusion is correct just for M > 0, as we are 
assuming. 
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FIG. 2. Left panel: /3 as a function of the event horizon for negative scalar charges. Depending on the value of the scalar 
charge, the number of solutions is one or three. Right panel; /3 as a function of the event horizon for positive scalar charge. 
Depending on the value of (3, there exist two or no solutions. Negative slopes correspond to thermally stable states. 



charge, we work out the temperature at the boundary, 
/3~^, as a function of r+ for fixed Q. This can be done 
by taking Eq. p^l): 



P{x,Q)^x- 



n/T^ 






where x 



-Itb. Q 



1 + (A-1) 
Q/r^ and ^ = 



(23) 



/3/ATrrjg. Observe 



that /3 positivity is guaranteed for Q > and is equivalent 
to require Eq. (151 for Q < 0. 

The plot of /3(a:, Q) is presented in Fig. |2]for specific 
values of Q. The curves are obtained for a certain A but 
the general features of the solution are independent of 
this choice. 

For the case Q > (right panel of Fig. I2| the function 
I3{x, Q) exhibits a simple behavior: if /? is below a cer- 
tain value /3maxj two black hole configurations with the 
same Q are allowed; otherwise no black hole (with the as- 
sumed scalar charge) is permitted. Moreover, the value 
of Z^max gets larger as Q decreases and it is bounded by 
the relation ^max < 2/'\/27 where the equality occurs for 
Schwarzschild black holes. 

For the case Q < (left panel of Fig. ^) /3 exhibits three 
different qualitative behaviors. For |Q| > |Qcrit| (with 
Qcrit < quantified later) at each temperature there is 
one single solution. In particular, the event horizon of 
the black hole grows as the temperature and reaches the 
boundary only when ^ = 0. For \Q\ < | Qcrit I there is an 
interval of temperatures outside which /3 behaves as in 
the previous case, i.e. only large (nearly extremal) black 
holes exist at very high (very low) temperature. Instead, 
for P within this interval, three black hole configurations 
are allowed. These three solutions approach the same 
event horizon for | Q| — > | Qcrit |~ and for this reason | Qcrit | 
can be defined as the value of Q at which the extrema of 
^, located at d/3/dx — 0, coincide. For instance, when 



A = 2 it turns out to be (see Appendix for details) 

1 



I Qcrit I = 



9-^4^5 



(24) 



B. The local stability 

The local stability in the saddle point approximation is 
related to the convergence of the integral in the on-shell 
partition function j26| , as we briefly review now. 

The on-shell partition function Z can be expressed as 



Z ^ f dE e-^^^^ ^ I dE v{E)t 



-0E 



(25) 



where i^(-E) is the density of classical states with en- 
ergy E. Applying the saddle point approximation, which 
consists in expanding the action up to quadratic or- 
der around the stationary points -Estat, the integration 
in Eq. (25 1 can be evaluated and the result is finite 
only if S^E^t^t > 0- On the other hand, in the sad- 
dle point approximation the energy and heat capacity 
can be approximated as {E) = —d\nZ/d(5 « Sstat and 
C = d{E)/dl3^^ sa dE^^^Jdp-^ and, finally, it turns out 
to be 



C 



P'{''AeJ 



(26) 



Consequently, since a configuration of a thermodynamic 
system is locally stable when its heat capacity is positive, 
one concludes that the convergence of the integral in the 
partition function is equivalent to the thermal stability 
of the system [26j . 

From this argument we can infer the stability property 
of the configurations depicted in Fig. [2] In fact, negative 
(positive) slopes of the isocharge curves /3(a;) are equiv- 
alent to positive (negative) heat capacity. Thus, for the 
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FIG. 3. The Euclidean action as a function of /3 for negative scalar charge in the three relevant cases. The values of Q axe the 
same as in Fig. [2] 



case Q > (right panel of Fig. I2| , when two black hole 
solutions are allowed at a given temperature, the smallest 
is thermodynamically unstable and the largest is stable. 
Instead, for Q < (left panel of Fig.j2| , when three black 
holes configurations with the same /3 and Q are possible, 
the smallest and largest are stable while the intermedi- 
ate is unstable. This unstable solution thus corresponds 
to a maximum of the Euclidean action whose value is 
associated to the tunnelling rate between the two stable 
configurations [37]. 

Therefore, the result presented in Fig. [2] might have 
striking implications for cosmology. By taking r-|_ <C rg 
we expect to mimic the black hole conditions in the 
present Universe, since the Universe is much larger than 
any event horizon radius. In such a case Fig. [2] shows 
that black holes with Q > have high temperature and 
tend to evaporate. Their relic abundance is then con- 
strained by the usual astrophysical bounds valid for ra- 
diating black holes [28] • On the contrary, black holes with 
Q < are cold long-lived configurations, which in general 
are poorly constrained and may be plausible dark matter 
candidates [2S]. Of course, the discovery of such stable 
objects would be an important finding in favor of massive 
gravity. Moreover, even though massive gravity can qual- 
itatively explain either dark energy and dark matter by 



tuning the gravitational force at long and intermediate 
distances _29j , the presence of black hole energy density 
might be a fundamental ingredient to fit the data. As 
a last remark, we highlight that for locally-stable black 
hole configurations, the action Ie fulfills the three re- 
quirements of Sect. In] Thus, for such configurations Z 
can be interpreted as the partition function of the system. 



From Eq. (22 1 one can therefore compute thermodynam- 



ical quantities as, for instance, the entropy: 






If 



1 , 



(27) 



which is the usual Bekenstein-Hawking formula [14] and 
agrees with the result found in Ref. [T3] by means of 
Wald's formula [30] ■ Hence, the Euclidean action re- 
duces to the free energy in the semiclassical approxima- 
tion, i.e. Ie = (iF with F = M - ThS, so that the first 
law of black hole thermodynamics is recovered by requir- 
ing F to have a minimum. Moreover, as demonstrated 
in Ref. [31] for static spherically symmetric gravitational 
systems, the thermal energy is identical to the quasilocal 
energy: 



£by = {E) = 



diE 

9/3 



(28) 
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FIG. 4. The Euclidean action as a function of /3 in the case Q = Q (left panel) and Q > (right panel). The solid and dashed 
lines correspond to a black hole and vacuum phase, respectively. 



Global stability, phase structure and critical 
behavior 



In the previous section we have determined when the 



black hole solutions ( 14 1 are locally stable. In particular, 



we found that under certain conditions multiple locally- 
stable solutions are allowed. However, for some fixed 
values of Q and /3, only one of these solutions can be 
globally stable (i.e. it corresponds to the global minimum 
of Ie) while the others have to be metastable (i.e. they 
are local but not global minima oi Ie)- 

In order to make manifest the phase structure of the 
black hole solutions (14 1, we analyze Ie as a function 
of temperature for a fixed Q. The result is shown in 
Fig. |3]and|4J respectively for negative and positive scalar 
charge. To read off the informations the figures contain, 
it can be useful to take into account also Fig. [2j When Q 
is negative and |Q| > |Qcrit| (top right and bottom pan- 
els of Fig. pi), there is only one solution independently 
of the specific temperature, as we already inferred from 
Fig. [2] Instead, when Q is negative but \Q\ < \Qcnt\ 
(top left panel of Fig. Isl), the competition between the 
black hole phases is more curious. Starting at very low 
temperature and then heating up the system we see that 
at the beginning only one solution exists. This phase re- 
mains globally stable up to the critical temperature /3~^ 
{(3c — 0.33 in the example of the Fig. ^ and wl) above 
which it becomes metastable. Subsequently, a first or- 
der phase transition may occur and, interestingly, above 
a certain temperature (/3 < 0.17 in the example) it has 
to. This property of the phase transition holds also start- 
ing at high temperature and then cooling down the sys- 
tem[^ Such a phase structure is isomorphic to the case of 



Reissner-Nordstrom black hole in AdS |33] and in a cav- 
ity [M] for the canonical ensemble (fixed electric charge) . 
To test the similarity, we examine in detail the behavior 
of our case at Q = Qcrit near the critical temperature 
Tc = /3~^. The entropy and the heat capacity scale as 
(see Appendix for details) 



S" - S'c ex (T - Tc 
C ex (T - Tc 



^-2/3 



(29) 
(30) 



implying that there is a second order phase transition at 
the critical temperature Tc- The value of the critical ex- 
ponent of the heat capacity has thus the same value as 
for Reissner-Nordstrom black holes in AdS, dS and flat 
space in a cavity at fixed electric charge [55H55] . Such a 
universal behavior is remarkable, knowing that the Eu- 
clidean action is different. 

In some sense, the phase structure for Q < and 
\Q\ < QcritI looks like the one for Q = (left panel 
of Fig. 4]) apart from two main differences: i) at high 
temperature the allowed phase is not a black hole solu- 
tion but the vacuum spacetime ([2| (dashed line in the 
figure); ii) starting at low temperatures and then heat- 
ing up the system, the first order phase transition (pos- 
sible at ;3 < 0.3 in the figure) is not always guaranteed 
since the vacuum phase is allowed at high temperaturerl 
Instead, for Q > (right panel of Fig. Ill the phase struc- 
ture is different: no phase is allowed at low temperature 
(/3 > ;5ijiax — 0.38 in the example). A possible expla- 
nation of this result might be that at low temperature 
the existing phase is a state that has scalar charge and 
no event horizon. This could perfectly be a star with 
low free energy [TTJ [H], which could not emerge from 



This feature might have interesting applications, as for instance 
in Randall-Sundruni or QCD models, which typically suffer from 
a too long-lived metastable phase |32| . 



^ Although the presence of the scalar fields, such a phase struc- 
ture is very similar to the one arising in GR for Schwarzschild 
black holes enclosed in a cavity |19| and Schwarzschild- AdS black 
holes [16]. 




0.4 0.6 0.8 1.0 





FIG. 5. Top left panel: /? as a function of the event horizon for a particular value of $ < $c- At a: = a;schw the scalar charge of 
the solution changes of sign. Top right panel: j3 as a function of the event horizon for a value of $ > $c. The black hole never 
engulfs the cavity and tends to the extremal limit at low temperatures. For x higher than the extremal limit, the configuration 
has no event horizon. Bottom panel: /3 as a function of the event horizon for different values of <l>. 



our analysis. On the other hand, one might guess that 
this pecuharity arises because we are assuming that Q is 
conserved. In the next section, we assume other bound- 
ary conditions where the scalar charge is allowed to vary 
inside the cavity. 



forward to see that the Euclidean action compatible with 



the regularity condition (171 is given by 



H£ 



BY 



4 



(31) 



where 



VI. OTHER BOUNDARY CONDITIONS 

In principle, one can assume the scalar charge not to 
be conserved inside the cavity. In such a case the phase 
evolution is no longer constrained by the initial value of 
Q but by the value of $, the scalar-charge "potential" at 
the surface of the cavity. In this ensemble it is straight- 



$ 



K 



,1-Al 



(32) 



One may prove that the present ensemble is always dom- 
inated by a configuration with a negative heat capacity 
and is therefore ill-defined, as we shortly see now. 

By combining Eqs. (18) and (32), /3 can be expressed 



as function of r_|_ and $. This relation can be more con- 
veniently rewritten as 



^ = 



X 



2+A(l_^A-l)2^ 



a;2(-l + A) - 2a;i+^(-l + A) - x^>'^^ + x^^ (-1 + A + $2 + (_i + a;)A$2) 

I 



(33) 



where x = r+/ri3, f3 = P/Atttb and $ = 2<I>rg ^. The behavior of the solution is plotted in Fig. [5] (solid lines) 
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for several values of <&. The scalar charge changes along 
the curve and it is positive (negative) when the solution is 
below (above) the dashed line. The radius at which this 
cross occurs is marked as Xgchw in Fig- p\ (top panels, the 
left plot being for $ < i>c and the right one for $ > $c)- 
Furthermore, when $ > $c, the black hole event horizon 
radius is bounded from above by the extremal case limit 
and can never reach the spherical cavity (see Appendix 
for some analytic results). 

The behavior of /3 described above seems problematic 
since the positiveness of the heat capacity depends on the 
sign of the slope of (3. For a positive slope, as in the case 
here, the heat capacity is negative. As we have discussed 
in Sect. VB this necessarily implies a divergent integral 



for the partition function, i.e. the partition function in 
the saddle point approximation is not well-defined and 
produces an imaginary result. This may be interpreted as 
the effect of an unidentified metastable phase. However, 
the present framework does not offer an unambiguous 
description and such interpretation should be taken with 
precaution. 



VII. CONCLUSIONS 

Lorentz breaking massive gravity is an interesting the- 
ory that may explain the recent acceleration of the Uni- 
verse without invoking dark energy. It also provides pe- 
culiar black hole solutions due to the presence of hair 
parameters that modifies the standard gravitational po- 
tential. In this theory the analog of the Schwarzschild 
black hole - the asymptotically flat spherically symmet- 
ric solution - depends on two parameters: the mass and 
the "scalar charge" (characterizing the hair strength). In 
this paper we have analyzed equilibrium states and phase 
structures of such a solution enclosed in a spherical sur- 
face kept at a fixed temperature. 

We have proven that when the scalar charge inside the 
cavity is not conserved, the ensemble is ill-defined. On 
the contrary, when the scalar charge is held fixed, the 
saddle point approximation can be applied to obtain the 
partition function Z = e~ '^ , where Ie is the regularized 
on-shell Euclidean action. With that formalism we were 
able to study the black hole thermodynamics and phase 
structure. In particular, the black hole entropy follows 
the Bekenstein-Hawking formula. 

Depending on the value of the scalar charge Q con- 
tained in the cavity, the phase structure presents com- 
pletely different behaviors that can be summarized as 
follows (being Qciit a critical value of the scalar charge): 

• Q > 0; Above a certain temperature there exist 
two black hole solutions with different event hori- 
zon radii. The smaller black hole is unstable while 
the larger is globally stable and describes the phase 
present in the cavity. Below that temperature we 
do not find any solution with event horizon, prob- 
ably because the analysis is not suitable for this 
phase. 



• Q — 0: Below a certain temperature no black hole 
solution exists and the phase consists of the usual 
globally-flat background of massive gravity. Above 
that temperature the previous phase competes with 
a new phase consisting of a large black hole. 

• Qcrit < Q < 0; There is an interval of temperatures 
where there are three black holes solutions with dif- 
ferent sizes of event horizon. The smallest and the 
largest are locally stable and therefore correspond 
to local minima of the free energy, but the interme- 
diate is unstable as it corresponds to a maximum 
of the free energy. Hence, a tunnelling between the 
two stable phases may occur and this flrst order 
phase transition is forced to happen when the tem- 
perature crosses the whole interval. 

• Q — Qciit-' Only one black hole solution is allowed 
at each temperature. At a certain temperature a 
second-order phase transition happens and at this 
moment the critical exponent of the heat capacity 
is -2/3. 

• Q < Qcrit-' There is a single (both globally and 
locally) stable black hole at every temperature. 

Further investigations on the cases with negative Q 
are worthwhile. Indeed, their isomorphism to Reisser- 
Nordstrom black holes in AdS, dS and flat space in the 
canonical ensemble |33H35j is peculiar, knowing that their 
Euclidean actions are different. Furthermore, at low tem- 
perature the globally stable configurations are small cold 
black holes. These objects overcome most of the as- 
trophysical constraints 528] and, in principle, might be 
the dark matter candidates of massive gravity. Dedi- 
cated dark matter analyses would be needed to check 
this possibility but, as a first step, one should under- 
stand whether the black hole scalar charge is actually 
conserved (or at least varies very slow) in nature. In 
order to address this question one should probably com- 
prehend the origin of the modified character of black hole 
horizons [101 136) , which is hard to understand without a 
known UV completion of the theory. However, with the 
use of the AdS/CFT correspondence one may attempt 
to study such problems |37j . In that context, the study 
of the phase structure may play an important role as 
a valuable test of the AdS/CFT correspondence. Since 
the phase structure emerging for AdS boundary condi- 
tions appears as well as for asymptotically-de Sitter black 
holes and asymptotically-fiat black holes in cavities |35) . 
one should expect the CFT side to have a rich variety 
of phases (such as deconfinement/confinement) dual to 
what has been studied in this paper. Moreover, Lorentz 
symmetry violation should arise in the CFT side |38) . 
Trying to understand regimes in which Lorentz violation 
is sizeable will hopefully allow us to devise the right ex- 
perimental tests to decide whether massive gravity is or 
is not realized in nature. 
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Appendix A: Some analytic results for A = 2 

In this appendix, we consider the case A = 2. This 
particular choice allows us to determine the main feature 
of the general case by some analytic calculations. In Ap- 
pendix |A 1| we focus on the ensemble with fixed scalar 
charge Q = Qcrit and in Appendix |A2| we determine the 
behavior of /3 for the ensemble with fixed scalar charge 
potential. 



The second derivative of /? vanishes, when evaluated at 
the critical point. From here, one can very easily obtain 
the behavior of the entropy and heat capacity near the 
critical temperature: 



S — Ocrit OC (T — Tcrit 
COC (T- Tcrit) 



-2/3 



(A7) 
(A8) 



As we have already mentioned in Sect. VC the critical 
exponent for the heat capacity has the same value —2/3 
as for Reissner-Nordstrom black holes in AdS, dS and flat 
space in a cavity at fixed electric charge [55H35) . This 
implies that at the critical value of the scalar charge, 
we have a second order phase transition. Moreover, this 
conclusion is independent of the value of A since the point 
a;crit is a stationary point of inflection as one can see from 
Fig-Il 



1. Case with fixed scalar charge Q 

For negative values of the scalar charge the function 
/9 has two extrema (see left panel of Fig. ^: one maxi- 
mum and one minimum. The location of such extrema is 
obtained by solving the equation 

5Q2 + (2 - 2Q - 3x)x^ - 6Qx{-l + Q + x)=0 . (Al) 

The critical scalar charge is obtained when both extrema 
are degenerate, i. e. w hen the discriminant of the poly- 
nomial equation ( Al ) is vanishing. By performing the 



explicit computation, the discriminant of the equation 
(Al) takes the form: 



A = -eoQ"E^"2" 



(A2) 



n— 3 



where S,i are real and positive constants. Therefore, sev- 
eral values for Q are allowed for a vanishing discriminant 
(A2). However, only one value of Q corresponds to a 



positive event horizon radius: 



1 



-9 - 4V5 



(A3) 



By replacing the critical value of the scalar charge in 
Eq. ( Al ), we obtain the critical value of the event horizon 



radius: 



■^crit 



= 5-2\/5 



By substituting both a;crit and Qcrit in P-. 
critical inverse temperature 



(A4) 
we have the 



38^5 . (A5) 

We note that such critical values are similar to the 
Reissner-Nordstrom case [3S]. At Q = Qcrit, we can ex- 
pand /3 around the critical point Xcrit^ 

1 ^3/3 



critical 



/3 - /3crit = 



3! dx^ 



(x - aicrit) + 



(A6) 



2. Case w^ith fixed scalar charge potential $ 

Depending on the value of the scalar charge potential, 
we have two qualitatively different behaviors for /3: 

• for < $ < l>c = 1: the extremal black hole limit 
is not allowed inside the cavity and the event hori- 
zon radius can therefore take all the possible values 
inside the cavity, i.e. < x < 1. At some temper- 
ature, the black hole engulfs the cavity. 

• for $ > i>c = 1: the black hole event horizon is 
bounded from above by the extremal black hole 
limit that takes place at a; = 1/$. Therefore, at 
very low temperatures the ensemble is dominated 
by nearly extremal black holes. In such a case, the 
black hole can never reach the cavity wall. 

In both situations, depending on the value of the event 
horizon radius, we can have positive or negative scalar 
charge configurations. The change of sign of Q takes 
place at 



'^schw 



1 
2I2 



l+4$2 
$4 



(A9) 



which is the value where the scalar charge vanishes, 
turning the hairy black hole into the conventional 
Schwarzschild black hole. Positive scalar charge config- 
urations are in the domain where x < Xgchw, implying 
therefore that black holes with negative scalar charge 
have X > Xschw 

At high temperatures, the ensemble is dominated by 
black holes with positive scalar charge. As the temper- 
ature drops, the scalar charge becomes smaller and the 
event horizon becomes larger. 
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